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We report on a joint experimental-theoretical study of collective diffusion in, and static shear 
viscosity of solutions of bovine serum albumin (BS A) proteins, focusing on the dependence on protein 
and salt concentration. Data obtained from dynamic light scattering and rheometric measurements 
are compared to theoretical calculations based on an analytically treatable spheroid model of BSA 
with isotropic screened Coulomb plus hard-sphere interactions. The only input to the dynamics 
calculations is the static structure factor obtained from a consistent theoretical fit to a concentration 
series of small-angle X-ray scattering (SAXS) data. This fit is based on an integral equation scheme 
that combines high accuracy with low computational cost. All experimentally probed dynamic 
and static properties are reproduced theoretically with an at least semi-quantitative accuracy. For 
lower protein concentration and low salinity, both theory and experiment show a maximum in 
the reduced viscosity, caused by the electrostatic repulsion of proteins. The validity range of a 
generalized Stokes-Einstein (GSE) relation connecting viscosity, collective diffusion coefficient, and 
osmotic compressibility, proposed by Kholodenko and Douglas [PRE, 1995, 51, 1081] is examined. 
Significant violation of the GSE relation is found, both in experimental data and in theoretical 
models, in semi-dilute systems at physiological salinity, and under low-salt conditions for arbitrary 
protein concentrations. 



I. INTRODUCTION 

A quantitative understanding of the dynamics in con- 
centrated solutions of interacting proteins is of impor- 
tance to the evaluation of cellular functions, and the im- 
provement of drug delivery. Transport properties such 
as collective and self-diffusion coefficients, and the static 
and high-frequency shear viscosities, are strongly affected 
by the aqueous environment [1], and in particular by 
crowding effects due to high concentration of macro- 
molecules, coupled both by direct and solvent-mediated, 
hydrodynamic interactions (His) [2-4] . The latter type of 
interaction, which is both long-ranged and of many-body 
nature, poses a particularly challenging task to a theo- 
retical treatment of diffusion and rheological transport 
properties. 

In the present paper, we report on a combined ex- 
perimental and theoretical study on collective diffusion, 
low shear-rate static viscosity, and static and dynamic 
scattering functions of concentrated solutions of bovine 
serum albumin (BSA) proteins. The goal of this study 
is twofold. On the one hand, we explore how far a sim- 
ple colloidal model in combination with state-of-the-art 
theoretical schemes can capture the microstructure and 
dynamics of proteins in solution. On the other hand, 
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we investigate the concentration- and salt-dependence 
of collective diffusion and the static shear viscosity, and 
use our results to test the validity range of a generalized 
Stokes-Einstein (GSE) relation which combines the col- 
lective diffusion coefficient with the isothermal osmotic 
compressibility and the shear viscosity. 

BSA is a protein which is readily soluble in water and 
stable over a wide range of salt and protein concentra- 
tions. Its stability and reproducibility make it well-suited 
as a model system of globular proteins. Proteins consti- 
tute identical solute units surpassing any synthetic col- 
loid suspension in terms of monodispcrsity. In this re- 
spect, they are ideally suited to the application of an- 
alytical theoretical models used with good success for 
large colloids. However, the construction of a quantita- 
tively accurate theoretical model for protein solutions is 
considerably obstructed not only by the potential pres- 
ence of impurities and oligomers, but also by the complex 
internal conformation and surface of a protein. The fold- 
ing state depends on various control parameters such as 
temperature, protein concentration, pH value, and salin- 
ity. The irregular protein surface implies an orientation- 
dependent protein interaction energy with repulsive and 
attractive parts, and furthermore complicates the de- 
scription of hydroynamically influenced transport prop- 
erties. 

In a first step towards calculating dynamic properties 
of proteins, it is nonetheless possible to use a model of 
reduced complexity, with system parameters such as the 
pH-dcpendcnt particle charge determined from a consis- 
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tent fit of theoretical expressions for the scattered in- 
tensity to the experimental static scattering functions. 
We use here a simple colloid model where the BSA in- 
teractions are described by the repulsive, electrostatic 
plus hard-core part of the isotropic Derjaguin-Landau- 
Verwey-Overbeek (DLVO) potential [5]. The effect of 
the non-spherical shape of BSA proteins is accounted for 
in the static intensity calculations within the so-called 
translational-oricntational decoupling approximation, by 
describing the proteins as oblate spheroids interacting by 
a spherically symmetric effective pair potential. 

Using this simplifying protein interaction model, the 
static structure factor, S(q), entering into the static scat- 
tered intensity, is calculated as a function of wavenum- 
bcr q, by using our newly developed modified penetrat- 
ing background corrected rescaled mean spherical ap- 
proximation (MPB-RMSA). This analytical method has 
been shown to be in excellent accord with numerically 
expensive computer simulation results for S(q) [6, 7]. 
The system parameters of the protein-interaction model, 
most notably the effective protein charge, are determined 
from adjusting the theoretically calculated static inten- 
sity, I(q), to the experimental one. The consistent agree- 
ment of calculated values and small-angle X-ray scat- 
tering (SAXS) data for I(q) in a wide range of concen- 
trations and wavenumbcrs indicates that left-out attrac- 
tive interaction contributions are of minor importance at 
the considered salinities. As an independent additional 
check, the static light scattering (SLS) data for S(q) at 
low q arc found to be well reproduced by the theoretical 
fits of the SAXS data. 

Without any further adjustment, the analytically cal- 
culated static structure factors are used as the only input 
to our theoretical calculations of the collective diffusion 
coefficient, dc, and the low shear-rate limiting static vis- 
cosity r\. To calculate dc and the high-frequency part, 
r)oo, of the static viscosity, we use two approximate an- 
alytical schemes, namely the pairwise additive hydrody- 
namic interaction (PA) approximation, and the so-called 
self-part corrected £7 method. As shown by two of the 
present authors [8] , these two methods give results which 
are in general in good agreement with more elaborate 
Stokesian Dynamics simulation results for particles with 
Yukawa- type pair interactions. 

The static viscosity, 

V = Voo + A77, (1) 

consists of a short-time part, 77^, determined solely by 
hydrodynamic interactions (His), and a shear-stress re- 
laxation part Ary, with A77 > 0. We calculate the latter 
using mode-coupling theory (MCT), which, like the two 
employed short-time schemes, requires S(q) as the only 
input. 

Our comparison with the experimental dc measured 
by dynamic light scattering (DLS), and with rj obtained 
from viscometry, is a stringent test for our theoretical 
results and for the employed isotropic interaction model, 
since except for the static input, no fit parameters are in- 



volved. In particular, no further adjustments of the the- 
oretical predictions have been made on referring to the 
actually non-spherical shape of BSA proteins. We show 
that despite the simplicity of our model, most dynamic 
features are well reproduced by the theoretical results, 
to an at least semi-quantitative accuracy. In particular, 
both a low-concentration maximum of the reduced vis- 
cosity, and a maximum in dc at a different concentration, 
are well captured by the theory. 

For BSA, also the short-time self-diffusion has been 
recently found to be reasonably well described by a sim- 
ple spheroid model [9, 10]. Of course, this does not im- 
ply that the complex conformation of a globular protein 
plays no role. The DLVO model (even with inclusion 
of van der Waals attraction) is not sufficient to fully ex- 
plain the rich phase behavior of proteins. For example, it 
has been shown that surface patchincss has an important 
effect on the phase diagram [11]. Also, binding of mul- 
tivalent ions to the protein surface can give rise to non 
mean- field behaviors beyond DLVO, such as charge in- 
version, re-entrant condensation and liquid-liquid phase 
separation [12-14]. 

Generalized Stokcs-Einstcin (GSE) relations, which 
approximately relate diffusion to rheological properties 
in concentrated complex liquids, are an important issue 
in microrheological studies, since a valid GSE relation 
allows to infer a rheological property more easily from 
diffusion measurements. Several GSE relations in col- 
loidal dispersions of electrically neutral (porous and non- 
porous) spheres, and charged particle suspensions have 
been explored [15-18]. We study here a GSE relation 
not discussed in this earlier work, which has been pro- 
posed by Kholodenko and Douglas [19]. This GSE rela- 
tion, which we refer to in the following as the KD-GSE 
relation, has been used in the biophysical and soft mat- 
ter community [20-23]. It relates dc to 77, and to the 
square-root of the isothermal osmotic compressibility. 

We present a thorough discussion of the validity range 
of the KD-GSE relation for BSA solutions, and for 
generic colloidal fluids of particles with screened Coulomb 
interactions, for a large range of salinities. Both the 
short-time and the long-time versions of the KD-GSE 
relation are considered. At high salinity, where the elec- 
trostatic interaction of particles is strongly screened, we 
find these two relations to become invalid at larger con- 
centrations. At lower salinity, the KD-GSE relations are 
poorly satisfied even at low concentrations. 

The paper is organized as follows. Sec. II includes the 
experimental details of the sample preparation, and of 
the SLS, DLS, SAXS, and rheological measurements. In 
Sec. Ill, we discuss the employed simplifying model of 
BSA, and present the essentials of our theoretical meth- 
ods, allowing for a fast calculation of measured static and 
dynamic properties. Our experiment data are shown in 
combination with the theoretical results in Sections IV 
and V, dealing with static and dynamic properties, re- 
spectively. Sec. V includes the examination of the KD- 
GSE relation. Our conclusions arc contained in Sec. VI. 
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II. EXPERIMENTAL DETAILS 



Static and dynamic light scattering 



Sample preparation 

BSA is a globular protein with a linear extension of 
about 7 nm. The considered aqueous solutions of BSA 
with no added salt, and with monovalent added salt such 
as NaCl, have a pH in between 5.5 and 7. Under these 
conditions, BSA is stable in solution, folded in its native 
state, and carrying a negative net charge in the range of 
roughly 8 to 20 elementary charge units (see below for 
details) [24, 25]. BSA was purchased from Sigma (cat. 
A3059) as a lyophilizcd powder, certified globulin- and 
protease free. 

The sample preparation for all experimental techniques 
started with the dissolution of protein powder in a sol- 
vent, and subsequent waiting until the solution was ho- 
mogenized. The protein mass concentration, c p , in the 
solution volume, Vn 2 o+iTip-9, is given by the BSA weight 



(2) 



Va 2 o 



with the specific protein volume 9 — 0.74 ml/g [2G] de- 
termines the self- volume of proteins upon dissolution. 

For small-angle X-ray scattering, dcionized and de- 
gased water was used as solvent. The samples with con- 
centrations higher than 15 mg/ml were prepared directly, 
while smaller concentrations were prepared from a stock 
solution of 18 mg/ml. The samples were filled into a 
plastic syringe and inserted into the capillary during the 
measurement. 

For the viscosity measurements, the solutions were pre- 
pared similarly using as solvent both deionized water, 
and solutions of NaCl in deionized water. The NaCl mo- 
larity is calculated from the total solution volume, in- 
cluding the protein self-volume. All solutions used for the 
viscosity experiments were further degased by a water-jet 
air-pump. 

For our light scattering experiments, stock solutions 
of BSA proteins in deionized water were mixed with so- 
lutions of NaCl in deionized water according to the re- 
quired concentration. The NaCl molarity is calculated 
from the total water volume. Then, every sample was 
pressed with a plastic syringe through a hydrophilizcd 
nylon membrane filter with a pore size of 100 nm (What- 
man Puradisc 13), and transferred into a cylindrical glass 
scattering cell. The cell was sealed immediately with a 
plastic cap. 

The effect of the difference in NaCl concentrations be- 
tween light scattering and viscosity samples, arising from 
the slightly differing sample preparation, is negligibly 
small. 



Multi-angle dynamic light scattering (DLS) was per- 
formed at various concentrations of protein and added 
salt, at a temperature of T = 295 K. In particular, the 
BSA mass concentration, c p , was chosen between 0.1 to 
150 mg/ ml, and the concentration of added salt was (no 
added salt), 5, 150 and 500 mM. Note that, even in the 
zero added-salt case, the analysis of the scattering data 
discussed in Sec. IV reveals a residual electrolyte concen- 
tration of a few mM, scaling roughly linear with c p (see 
Table I). This suggests a few possible sources of the resid- 
ual electrolyte ions. Firstly, a possible source could be the 
surface-released counterfoils of charged BSA oligomers, 
not contained in our monodisperse model. Secondly, a 
salt contamination of the BSA stock, and thirdly the dis- 
sociation of acidic or alkaline surface groups off the BSA 
proteins cannot be excluded. 

Static light scattering (SLS) experiments were per- 
formed on the same samples. We used a combined 
SLS/DLS device from ALV (goniometer: CGS3, correla- 
tor: 7004/FAST), located at the Institut Laue Langevin 
in Grenoble, with a minimum correlation time of 3.125 ns 
as initial and shortest time. The HeNe laser was operat- 
ing at wavelength Ao = 632.8 nm, with an output power 
of 22 mW. The accessible range for the scattering angle 
(wavenumber) was 30 - 150° (q = 0.007 - 0.026 nm" 1 ). 
Moreover, the DLS intensity autocorrelation function de- 
cays on a time scale much slower than the interaction 
time, rj ~ cr^ /(4do) ~ 0.3 /Lts, of BSA, where do is 
the single protein average translational free-diffusion co- 
efficient, and (J do is an effective hydrodynamic diameter. 
Hence, DLS probes the long-time collective diffusion of 
BSA, in the q ->• limit. 

The normalized intensity autocorrelation function ob- 
tained from DLS, 



(Z(q,0)I(q,t)) 

</(q)> 2 



was fitted according to the Siegert relation, by the double 
exponential decay function 



g 2 (q,t)-l= f A--exp[-A<7 2 *]J +B, (3) 

with decay constants D\ and D2, and amplitudes A\ and 
A-}. The fit results were essentially the same with and 
without the background-correction constant B. At all 
probed angles, the two decay constants are widely sepa- 
rated (Di 3> D2). The faster mode, D±, is attributed to 
the (long-time) collective diffusion coefficient, dc, of BSA 
monomers. The appearance of the slower mode charac- 
terized by L>2, can be attributed to the slow motion of 
the larger impurities and oligomers. After having checked 
that D\ is overall g-independent within the experimen- 
tal resolution, it was averaged with respect to its resid- 
ual scattering angle fluctuations to gain better statistics. 
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Data on Di are rather noisy in comparison to D\, and 
show no clear dependence on q, c p , and on the concen- 
tration of added salt. 



Small-angle X-ray scattering 

Aqueous solutions of BSA with mass concentrations 
between 0.9 mg/ml and 270 mg/ml, and without added 
salt, were measured by small-angle X-ray scattering 
(SAXS), at the beam line ID02 of the European Syn- 
chrotron Radiation Facility (ESRF) in Grenoble, France. 
The standard configuration at a 2m samplc-to-dctector 
distance, and a photon energy of 16051 eV was used. 
Measurements were repeated several times in the flow 
mode and with short detection times to ensure the ab- 
sence of radiation damage. The data from the CCD 
were processed with the standard routines available at 
the beam line for radially averaging the data and cor- 
recting for transmission. Repeated measurements were 
summed up, and the solvent scattering was measured in- 
dependently and subtracted from the data. Additionally, 
two dilute samples (at c p = 1 and 2 mg/ml) with 150 mM 
of added NaCl were measured for form factor fitting. 



Viscosity measurements 

The viscosity data were measured at T = 25° C, for 
different concentrations of protein and added salt. The 
first dataset was obtained for solutions without added 
salt, while the second set describes systems with 150 mM 
NaCl. All measurements were performed at a shear rate 
of 60 Hz <C I/17, using the suspended couette-type vis- 
cometer described in Ref. [27]. The important advan- 
tage of this instrument is the possibility to collect data 
without errors caused by the surface shear-viscosity. A 
test made for c p « 20 mg/ml and 100 mg/ml, without 
salt and for 2M added NaCl, revealed no shear-rate de- 
pendence of the viscosity for shear rates between 50 and 
95 Hz. The precision of the viscosity measurements is 
approximately 0.1%. In order to minimize systematic 
errors, every measurement was repeated three times, in- 
cluding separate sample preparations. 

The viscometer directly measures the relative shear- 
viscosity of the solution against pure water (for technical 
details see Ref. [27]). For the aqueous BSA solutions 
without added salt discussed in this work, the relative 
viscosity was directly measured. For BSA solutions with 
added salt, this quantity was obtained as the ratio of the 
following two values: (a) the directly measured relative 
viscosity of the BSA solution with salt against water di- 
vided by (b) the directly measured relative viscosity of 
the salt solution (without BSA) against water. 



III. THEORY 
Single-particle properties 

In the following, we discuss the spheroid model of BSA. 
We use this model for the form factor fitting, and in de- 
termining effective sphere diameters related to different 
single-particle properties. 

At low protein concentration and sufficient amount of 
added salt, inter-protein correlations arc negligible. The 
scattered intensity, I(q), is then solely determined by the 
form factor P, i.e. /(g) oc P(q). Crystallographic mea- 
surements [28-30] have revealed a flat and roughly heart- 
shaped structure of albumins. The computation of single- 
particle properties with an account of the highly complex 
particle shape of biomolecules can be done by numerical 
simulations only and is beyond the scope of this paper 
[29, 31]. Rather, the aim of the present study is to give 
an essentially analytic description of the microstructure 
and the dynamics of interacting BSA proteins with low 
computational cost. We therefore intentionally choose an 
extremely simple model for the fit of the protein form fac- 
tor, by an oblate, solid ellipsoid (spheroid). Clearly, this 
mapping of the complex protein configuration onto an 
essential geometric shape is a delicate and broad topic 
on its own. Considering that the focus of the present 
work is on collective correlations rather than on single- 
particle properties, we cannot discuss all details of this 
subtle matter; we basically follow the approach of Ref. 
[32]. 

For a homogeneously scattering spheroid with dimen- 
sions a and b, where a denotes the semi-axis of revolution, 
the orientationally averaged form factor, P e u, is given by 
[33] 

PeiM) = [ dn\f{q,n)\ 2 (4) 
Jo 

with the scattering amplitude f(q,n) = 2>j\{u)/u, and 
u = q^Ja 2 y? + b 2 (l — fi 2 ). Here, j\ is the spherical 
Bessel function of the first kind. 

The fit of Eq. (4) to our newly recorded, low- 
concentration SAXS intensities at c p = 1 and 2 mg/ml, 
and for 150 mM of added NaCl, is shown in Fig. 1 of 
Sec. IV, along with a discussion of the obtained best fit 
values a = 1.75 nm and b = 4.74 nm. 

When protein correlations come into play at higher 
concentrations or lower salinities, the spheroid model of 
BSA becomes too complex for an analytic treatment. 
Therefore, as far as the protein-protein interactions are 
concerned, we describe the proteins as effective spheres 
with diameter a. Depending on the considered single- 
particle property, different definitions for a can be given. 

Consider first the geometric effective diameter, a geo = 
8(a6 2 ) 1 / 3 = 6.80 nm, which follows from equating the 
volume of the effective sphere to that of the spheroid. 
This effective diameter reflects the volume of the pro- 
tein and the hydration layer visible to SAXS, but does 



5 



not include thermo- and hydro-dynamic effects of non- 
sphericity [34, 35]. Thus, it should be considered as a 
lower boundary to the effective sphere diameter. 

A thermodynamic effective diameter, as 2 = 7.40 nm, 
follows from demanding equal second virial coefficients, 
£>2{T), of hard spheroid and effective hard sphere [36]. 

Alternatively, dynamic single-particle properties can 
be used in defining the effective diameter. For hydro- 
dynamic stick-boundary conditions and a < b, the trans- 
lational free diffusion coefficient of an isolated spheroid 
reads [29, 37, 38] 



df{a,b) 



ksTS(a, b) 
12irr]oa 



(5) 



with absolute temperature T, Boltzmann's constant ks, 
solvent shear-viscosity rjo, S(a, b) = 2 atan £(a, &)/£(a, b) 
and £(a, b) = y/\a 2 — b 2 | /a. Equating dff 1 to the dif- 
fusion coefficient, c?o = /cBT/(37r?7o<7d ), of an effective 
sphere gives Od a = 7.38 nm. 

Finally, one can derive another effective diameter from 
the intrinsic viscosity 



[, ? ]=hm!M_^ 
0->o 7700 



(0) 



where (f> is the particle volume fraction. For a spheroid 
with hydrodynamic stick-boundary conditions [39, 40], 



{] ell = 1 + ^L 

° 15tt 
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- - f 


- 0.628 




a 




1 



1-a/b 



0.075a/6 



(7) 



which for a = b reduces to the Einstein result, [rj\ sp = 
2.5, for a solid sphere. Note here that [rff 11 > 2.5 for a ^ 
b. Explicitly, [rff 11 = 3.25 for the best fit values a and b 
given in Fig. 1. On demanding equality of the interaction- 
independent linear terms in the virial expansions of the 
viscosity, 

^- = 1 + [t#+O(0 2 ), 
Vo 

for spheroids and effective spheres, and on using cj) eU = 
(47r/3)a6 2 n and <fi 3ph = (n /fya^n for an equal num- 
ber density n, the effective diameter crui = 7.42 nm is 
obtained. 

Since the aspect ratio, b/a = 2.71, is rather close to 
unity, the four obtained effective diameters are quite sim- 
ilar in magnitude. We use a = ub 2 =7.40 nm in all our 
calculations of static and dynamic properties discussed 
in this paper. 



Static scattering intensity and structure factor 

Concentrated protein solutions exhibit pronounced 
inter-particle correlations which are reflected in the static 
scattering intensity. This applies also to dilute, low- 
salinity solutions where the proteins show long-ranged 
electrostatic repulsion. 



In order to allow for an analytical theoretical treat- 
ment, we assume that the static scattering intensity of 
interacting BSA proteins can be approximated by 



I(q) = AcpP e a(q)S m {q), 



(8) 



where S m is the so-called measurable static structure fac- 
tor. Here, A is a q-independent factor (of dimension 
velocity 3 ), that should be the same for all intensity mea- 
surements corrected for recording time and source inten- 
sity. 

For calculating S m (q), we use the rotational- 
translational decoupling approximation [41, 42], where 
the spheroid shape is accounted for in the scattering am- 
plitudes only, so that 



S m (q)= l-X(q) +X(q)S(q) 



Here, 



X(q) 



Pell(q) 



d^f(q^) 



(9) 



(10) 



with < X(q) < 1 and X(q -> 0) = 1, and S is the so- 
called ideal structure factor of ideally monodisperse effec- 
tive spheres of diameter a = ub 2 and screened Coulomb 
repulsion of DLVO type. For the BSA model spheroid 
used here, X(q) stays close to unity for q < 0.5 nm -1 , 
decaying for larger q steeply towards its first zero value 
atgw 1.3 nm -1 . For q > 1.3 nm" 1 , X(q) < 0.04. The 
orientational disorder assumed in the decoupling approx- 
imation has the general effect of damping the oscillations 
in S m (q). While S m (q) is practically equal to one for 
q > 1.3 nm -1 , irrespective of the still visible oscillations 
in S(q), the effect of orientational disorder on S m (q) is 
weak in the range q < 0.5 nm -1 , where the most dis- 
tinctive features in S(q) arc seen. We further note that 
S m (q — > 0) = S(q — > 0) for monodisperse systems, a 
feature which plays an important role in our upcoming 
discussion of collective diffusion. 

The ideal structure factor, S(q), entering into Eq. (9), 
is calculated using the repulsive part of the DLVO pair- 
potential [5], 

{oo , x = r/cr < 1, 
e~ kx (11) 
7 , x>l, 
x 

also referred to as the hard-sphere Yukawa (HSY) po- 
tential. The coupling parameter, 7, and the screening 
parameter, fc, are given by 

k 2 = k\ + kl = (2A<t>\Z\ + 8irn s a 3 ) . (12b) 

1 — 

Here, Lb — fie 2 je is the solvent-characteristic Bjerrum 
length in Gaussian units, (3 = l/(/csT), e is the solvent 
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dielectric constant, and Z is the effective protein charge 
number in units of the proton elementary charge e. The 
factor 1/(1 — (f>) ink 2 corrects for the free volume available 
to the microions [43, 44]. We have not included van der 
Waals (vdW) forces in u(x). However, we have checked 
that the influence of vdW attractions is small for most 
of the considered systems. 

Eq. (12b) consists of two additive parts. The first part, 
k 2 oc \Z\, is due to protein-surface released counterions, 
which are assumed to be monovalent. The second part, 
k 2 , accounts for the screening due to all other monova- 
lent microions. Owing to the overall charge neutrality, 
this contribution is proportional to the co-ion concentra- 
tion n s . A lower bound of n s > 10~ 7 M in pH-neutral 
aqueous solutions is due to the self-dissociation of water. 
Additional contributions to n s can arise from dissolved 
CO2, and added salt such as NaCl. For a protein so- 
lution, n s can have a (putatively linear) dependence on 
c p if charged protein oligomers are present, acting as an 
additional source of surface-released counterions not con- 
tained in our model. Moreover, the protein stock solution 
might contain a residual amount of salt, and the proteins 
might dissociate acidic or alkaline surface groups during 
solvation. Note that due to the overall charge neutral- 
ity, the total concentration of monovalent counterions is 
given by n s + 6(j)\Z\/ (na 3 ). 

In recent work [G, 7], two of the present authors 
have derived a computationally efficient integral equation 
scheme for computing S(q) using the screened Coulomb 
potential in Eq. (11). This so-called modified penetrating 
background corrected rescaled mean spherical approxi- 
mation (MPB-RMSA) shares the analytical simplicity of 
the widely used RMSA [32, 45], but is distinctly more 
accurate. All calculations of S(q) in this paper are based 
on the MPB-RMSA. 

The spheroid- Yukawa (SY) model used in our cal- 
culations of I(q) and S m (q) ignores orientational- 
translational coupling. Therefore, it can be expected 
to apply only to fluid-phase BSA solutions when c p is 
sufficiently low, and when the ionic strength is not too 
large, so that the anisotropic protein shape and pair- 
interaction parts are not important. At larger c p , there 
is orientational-translational coupling, and the decou- 
pling approximation becomes invalid. We note again 
that the possible presence of residual protein oligomers 
and scattering impurities is not accounted for in our one- 
component model. The virtue of the SY model, however, 
is its analytical simplicity. The concentration range in 
which the SY model is applicable to BSA is examined in 
Sec. IV. 

Since we use a spherically symmetric screened 
Coulomb plus hard-core pair potential for the protein- 
protein interactions, a short discussion of the neglected 
anisotropy in the electric double layer around a charged 
spheroid is in order here. 

The mean electrostatic potential, <I>(r, /1) = 
$;(r)P;(/i), of a spheroid with a correspond- 
ing axisymmetric charge distribution immersed in an 



electrolyte solution includes in general higher-order 
multipoles with / > 0. Here, r is the distance of the 
spheroid center to the field point, [i = cosi? is the cosine 
of the angle relative to the spheroid rotational symmetry 
axis, and the P/'s are Legendre polynomials. 

For large r, all multipoles decay asymptotically equally 
fast according to [46-51] 

p -rer 

Mr) ~ fi , (13) 

r 

where k denotes the inverse electrostatic screening 
length, and /; depends on the charge distribution. This 
implies that, in principle, the pair interaction energy of 
two spheroids depends on their relative orientation even 
when r 3> However, the multipolar strengths, 

for a spheroid with b/a ~ 1 can be expected to be small 
for larger I. Moreover, since after orientational averag- 
ing, (Pi(/j,)) = for all I > 0, our neglect of anisotropic 
pair interaction contributions can be expected to be rea- 
sonable, for systems where the particles can essentially 
rotate freely. 

Short-time diffusion 

We summarize here the analytical methods used in cal- 
culating the (short-time) collective diffusion coefficient 
dc- These methods require S(q) as their only input, with 
the BSA protein interactions described by the spherical 
pair potential in Eq. (11). 

The colloidal short-time regime covers correlation 
times t within rg C t < rj. Here, tb = m p /(37rr)oa) 
is the momentum relaxation time of a globular protein 
of mass m p . Within a short-time span, a protein has dif- 
fused a very small fraction of its size only. For BSA in 
water, tb ~ 1 ps, and 77 ~ 0.3 fis. The BSA short-time 
dynamics is thus not resolved in our DLS experiment 
determining the measurable dynamic structure factor, 
S m {q,t), as a function of wavenumber q and correlation 
time t. 

Within the translational-oricntational decoupling ap- 
proximation used in the SY model, S m (q, t) is determined 
by the right-hand-side of Eq. (9) with S(q) replaced by 
S(q,t). The latter is the ideal dynamic structure factor 
of ideally monodisperse, charged effective spheres inter- 
acting according to Eq. (11). 

Owing to the smallness of the proteins compared to 
the wavelength of visible laser light used in our DLS ex- 
periments, one obtains f > t; and q <C q m - Here, q m is 
the wavenumber where S attains its principle peak value. 
Since X(q <C q m ) ~ lj it follows that S m (q <C q m ,t) ~ 
S(q <C q m , t), so that the influence of orientational disor- 
der on the measured S m (q, t) via the spheroid form factor 
is negligible. 

As a consequence, DLS determines the long-time col- 
lective diffusion coefficient, d^, according to 

S m (q < q m ,t > 77) oc exp [~q 2 d^t] . (14) 
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The coefficient d c , also referred to as the gradient diffu- 
sion coefficient, quantifies the long-time decay of long- 
wavelength, isothermal protein concentration fluctua- 
tions. In Eq. (14), additional scattering contributions 
to S m (q,t), originating from oligomers and large impu- 
rities, are neglected. As discussed in relation to Eq. (3), 
these give rise to an additional, exponentially decaying 
mode with a mean diffusion constant, which is sub- 
stantially smaller than d c . 

While, in principle, d c needs to be distinguished from 
its short-time counterpart d c , with d c < d c , it has been 
shown [41, 52] that the relative difference is very small (< 
5%) even in highly concentrated systems. For solutions 
like the ones considered in this work, where non-pairwise 
additive HI contributions are small, dc = d c becomes 
practically identical to d c . This allows us to use more 
simple short-time dynamic methods for calculating dc- 

To this end, we use two complementary analytical 
methods, namely a self-part corrected version of the so- 
called 8 j scheme due to Bccnakkcr and Mazur [8, 53-56], 
denoted here as the corrected £7 scheme for brevity, and 
a pairwise additive (PA) approximation of the His. The 
latter becomes exact at very low concentrations, but its 
prediction for d c worsens when protein volume fractions 
<f> > 0.05 are considered (see our discussion of Fig. 3 
in Sec. V). On the other hand, the PA predictions for 
7700 , and for the short-time self-diffusion coefficient d$ not 
considered here, are reliable up to substantially larger 
volume fractions, as has been ascertained in comparison 
to Stokesian Dynamics computer simulations [8, 17] and 
experimental data [53]. The PA expression for d c reads 

/oo 
dxx 2 g(x)yf 2 (x) 

+ 80^ dxx 2 g(x)[x a 12 (x) -j/? 2 (:c)]}, (15) 

with ds given in PA approximation by 



d 



dxx 2 g(x)[xi 1 (x) + 2yt 1 (x)-3]. (16) 



The two-body mobility functions, x^ and y^-, can be ex- 
panded analytically in powers of a/r = 1/x. The short- 
range mobility parts 



Xi 2 (x) 

y? 2 (*) 



Xi 2 (x) 



(3/4) af 
(3/8) aT 



(1/8) x- 3 , 
(1/16) x~ 3 , 



include all terms in the series expansion in 1/x with the 
far-field terms up to the dipolar level subtracted off. For 
x > 3, an explicit analytical expansion to 0(x~ 20 ) is 
used [57]. Since the series expansion in 1/x converges 
slowly at small separations, accurate numerical tables, 
which account for lubrication at near-contact distances 
[58], are employed for x < 3. 



The only input required in Eqs. (15) and (16) is the 
radial distribution function g, related to S by a one- 
dimensional Fourier transform [59]. The two functions 
are obtained in our analysis by the analytical MPB- 
RMSA. 

The second short-time method used in the present 
work for calculating d c « dc and 7700, is the self-part 
corrected £7 scheme. In this scheme, d c is obtained from 
the exact relation [41] 



4^ 



, lim H d (q) 
do ?->o 



(17) 



containing the distinct part, Hd(q), of the so-called hy- 
drodynamic function H(q). The ^7 scheme of Bccnakkcr 
and Mazur provides an easy-to-use integral expression for 
Hd(q), including S(q) as the only required input. The ex- 
plicit form of the 57-scheme expression for Hd{q) is given 
in [17, 56] and will be thus not repeated here. 

Extensive comparisons with Stokesian Dynamics simu- 
lations [17], and experiments on charged colloids [53, 60], 
and for small (f> also with PA calculations, have shown 
that the 67 scheme predictions for Hd{q) arc quite good 
for all concentrations up to the freezing transition value, 
even though the £7 scheme involves hydrodynamic ap- 
proximations at any concentration. In particular, it dis- 
regards lubrication effects. Lubrication, however, is in- 
consequential for charge-stabilized particles where near- 
contact configurations are unlikely. 

Different from Hd(q), the accuracy of the #7 scheme 
is less good for charged particles regarding the self-part, 
ds, of d c in Eq. (17) [17, 53]. To remedy this deficiency, 
we use a hybrid method, referred to as the self-part cor- 
rected scheme, in which ds is calculated using the 
PA expression in Eq. (16). It has been shown both for 
charged colloids [8, 17, 53] and Apoferritin protein solu- 
tions [til], that this hybrid method works quite well at 
fluid state concentrations. 



High-frequency viscosity 

The high-frequency viscosity, 77^, linearly relates the 
average suspension shear stress to the average rate 
of strain in a low-amplitude, high-frequency oscilla- 
tory shear experiment. While this short-time quantity 
has been rather routinely determined for micron-sized 
charge-stabilized colloids [16, 62], a direct mechanical 
measurement of for BSA solutions is difficult, since 
the required frequencies to 3> tJ x are in the MHz regime. 
We are interested here in rjoo since, according to Eq. (1), 
it is an important contribution to the static viscosity 77. 
The latter has been determined experimentally in the- 
present work. 

In PA approximation, 77^ is given by [17, 63, 64] 

^ = 1 + -0(1 + 0) + 6O0 2 / dxx 2 g(x)J(x), (18) 
Vo 2 J 1 
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where the rapidly decaying shear mobility function J(x), 
with J(x) = 15/128 x- 6 + 0(x" 8 ) for stick boundary 
conditions, accounts for two-body HI effects. In per- 
forming the integral over g(x), the leading-order long- 
distance contribution is dominating for x > 3. Accurate, 
numerical tables, where the lubrication effect for x ~ 1 
is included, are used for x < 3 [58]. 

The Sj scheme of Beenakker and Mazur can be also 
used for calculating r/oc. Similar to the ^-scheme ex- 
pression for d^, the standard (2 nd order) £7 scheme re- 
sult for 7700 consists of a microstructure-independent self- 
part, and a distinct part given in form of an integral over 
S(q) [55]. In recent work, two of the present authors have 
shown that a self-part corrected version of the original $7 
scheme expression for gives results for charged par- 
ticles in very good agreement with Stokesian Dynamics 
simulations [8]. This self-part corrected £7 scheme for 
rjoc is used in the present work. 



Static shear-viscosity 

In long-time rheological measurements on protein so- 
lutions under steady shear, there is an additional shear- 
stress relaxation part, A77, contributing to the static vis- 
cosity r) = r/oo + Ary. This contribution is influenced 
both by His and direct interaction forces. It can be cal- 
culated approximately within the mode-coupling theory 
(MCT) of Brownian systems. While a version of MCT 
for Arj with far-held HI included has been discussed in 
earlier work together with an extension to multicompo- 
nent systems [65], for analytical simplicity we use here 
the standard one-component expression 

(19) 

which has been obtained, e.g. in [65], under the ne- 
glect of His. In principle, Arj MCT should be calculated 
self-consistently by a numerically expensive algorithm in 
combination with the corresponding MCT memory equa- 
tion for S(q, t) [66]. However, the BSA solutions explored 
here are rather weakly coupled particle systems, with 
structure factor maxima S(q m ) < 1.2. Thus, as we have 
thoroughly checked in comparison to fully self-consistent 
MCT calculations, Ar/ MCT can be obtained more sim- 
ply in a first iteration step where S(q,t) in the inte- 
gral of Eq. (19) is approximated by its short-time form 
S(q,t)/S(q) = exp[-q 2 d t/S(q)}, valid without HI. The 
difference to the fully self-consistent result for Aij MCT is 
at most a few percent, even for the most concentrated 
systems considered. 

Moreover, again due to the only moderately strong in- 
terparticle correlations, Arj augments 77^ by at most ten 
percent. Therefore, the neglect of HI in Arj MCT can be 
expected to be rather insignificant for the systems consid- 
ered since the dominant effect of HI is included already 
in r/oo. Theoretical results for r/ shown in this paper are 




Fig. 1. BSA form factor fit. Open circles: SAXS intensities at 
two protein concentrations of c p = 1 and 2 g/l, for 150 mM of 
added NaCl. The SAXS intensities have been divided by c p , 
and by a common, ^-independent factor Aq. Red solid line: 
Angular-averaged spheroid form factor according to Eq. (4), 
fitted to the SAXS data within 0.3 nm -1 < q < 1.15 nm - , 
as indicated by the blue vertical line segments. The obtained 
fit values are a — 1.75 nm and b — 4.74 nm. Inset: Intensity 
on a double linear scale. 

all based on the first iteration solution for Arj MCT , and 
on 7700 calculated using the self-part corrected £7 or PA 
schemes. For all explored systems, the difference in rjoo 
between the PA and corrected 5j scheme is at most two 
percent. 

IV. STATIC PROPERTIES: EXPERIMENT AND 
THEORY 

A. Form factor fit 

In Fig. 1, SAXS intensities for BSA solutions of very 
small protein weight concentrations, c p = 1 and 2 mg/ml, 
and 150 mM of added NaCl, are shown along with the 
best-fit spheroid form factor. Note that our form fac- 
tor fit relies on a simplified shape model, so that some 
controlled systematic deviations from experimental data 
are to be expected. To check for a residual effect of in- 
terparticle correlations on I(q), S(q) was calculated for 
the present two systems to first order in </> using the full 
DLVO potential, with |Z| ~ 30 and a Hamaker constant 
of 3/cbT [67]. The so-obtained structure factor deviates 
only very little from unity with S(q —> 0) « 0.99. If vdW 
attraction is ignored, S(q -> 0) is slightly lowered to 0.98. 

Thus, to fit the measured intensity in Fig. 1, we have 
used Eq. (8) for I(q) with S m (q) set equal to one. Us- 
ing an automatic weighted least-squares minimizer, the 
spheroid semi-axes a and b entering into P e u (q) were var- 
ied to achieve a best fit intensity for a given prefactor A 
in Eq. (8). This fitting procedure was iterated for differ- 
ent values for A, until optimal agreement with the SAXS 
intensities within the range 0.3 nm -1 < q < 1.15 nm -1 
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was achieved, resulting in a = 1.75 nm and b = 4.74 
nm. These values for the spheroid semi-axes are in good 
accord with previously reported values, and in reason- 
able agreement with the linear dimensions of the reported 
heart-shape like crystal structure of albumins [28-30, 32]. 
In a related, recent study by part of the present authors 
[9], similar values a = 1.80 ±0.05 nm and b = 4.60 ±0.15 
nm have been determined, which are in decent agreement 
with the values obtained here. The optimized value for 
A, denoted by A , has been also used in our SAXS in- 
tensity fits for systems without added salt, which will be 
discussed in the following subsection. 

The best-fit form factor, P e u, depicted in Fig. 1 de- 
viates from the SAXS intensities outside the fitted g- 
range. For q > 1.15 nm -1 , corresponding to length scales 
2ir/q < 6 nm < <j, the complex internal structure of BSA 
is probed, which is not accounted for in our simplify- 
ing SY model. The deviations visible for q < 0.3 nm -1 , 
corresponding to distances of roughly 20 nm or larger, 
are likely due to additional scattering species made up 
of larger particles such as BSA oligomers or impuri- 
ties. Since the size-, form-, and charge-distributions of 
oligomers and impurities are unknown, our choice of the 
lower q-boundary in fitting /(g) is somewhat more am- 
biguous than the upper boundary. Therefore, we have 
repeated the intensity fitting for various low-q bound- 
aries, finding that the weighted least squares deviation 
increases dramatically if the boundary is selected below 
0.3 nm -1 . Moreover, the fit values for a and b remain 
essentially constant when the lower o-boundary is chosen 
larger than 0.3 nm" 1 . 

The fit parameters of a spheroid form factor to SAXS 
data of proteins in general depend slightly on the mea- 
sured q range, the prepared protein concentration, sol- 
vent and salt conditions, and background subtraction. 
In the context of the present study, the related changes 
of the spheroid model parameters are small compared to 
the experimental error bars and will be discussed in the 
next section. 



B. Concentration series of scattered intensities 

Fig. 2 includes the SAXS intensities for all explored 
BSA solutions without added salt that could be fitted us- 
ing the decoupling approximation expression in Eq. (8), 
for S(q) calculated in MPB-RMSA using the screened 
Coulomb potential in Eq. (11). In order to emphasize 
the shape differences across the dataset, the intensities 
are divided by their respective fitted amplitudes A, and 
by the protein concentrations c p . The most concentrated 
solution shown here is the one for c p = 90 mg/ml. Two 
even more concentrated systems for c p = 180 and 270 
mg/ml are not depicted in the figure since their intensi- 
ties could not be fitted reasonably well by the SY model. 

In order to fit the experimental intensity data using 
Eq. (8) , some deviations of the prefactor A from the op- 
timized form factor fit value Aq have to be allowed for 




0.3 0.6 ] 0.9 1.2 

q [nm ] 



Fig. 2. SAXS intensities from BSA solutions at various Cp, 
without added salt, divided by Ac v . From top to bottom: 
c p = 0.9,1.8,4.5,7.2,9,13.5,18,45, and 90 mg/ml. The in- 
tensity curves are displaced in steps of 0.1 along the vertical 
axis for better visibility. The SAXS-data for the extended 
range q < 4 nm -1 were taken into account in all fits, but 
shown here only up to q — 1.2 nm - . Red solid lines: best 
fits according to Eq. (8) with S(q) calculated in MPB-RMSA. 
The fit parameters are listed in Table I. 



(see Table I). The fit of each individual intensity curve 
in Fig. 2 was made as follows: After dividing the SAXS 
intensity by Aq and c p , the weighted sum of quadratic 
deviations between SAXS data points and the intensity 
according to Eq. (8) was minimized by an automatic 
three-dimensional weighted least-squares minimizer with 
respect to the fitting parameters {\Z\,n s ,<fi}. For each 
concentration, the whole experimental dataset was used, 
for wavenumbers from 0.07 to about 4 run . If the fit 
was unsatisfactory, the prefactor A was slightly altered, 
and the optimization with respect to {\Z\, n s , </>} was re- 
peated. This procedure was iterated until convergence 
in all fit parameters was achieved. For all considered 
concentrations, Lb = 0.711 nm, a = (tb 2 = 7.40 nm, 
a = 1.75 nm, and b = 4.74 nm were kept fixed. Table I 
summarizes the obtained best fit parameters. 

While the overall intensity fits for the two lowest con- 
centrations, c p = 0.9 and 1.8 mg/ml, look quite reason- 
ably good, they contain some peculiarities. A shoulder 
is present in the fit intensity extending from q rj 0.3 to 
0.8 nm -1 , overshooting the experimental data by several 
standard deviations. Moreover, the prefactor A is sub- 
stantially larger than A in both cases, and the fitted 
effective charge number \Z\ assumes a questionably large 
value of 34.5 for c p = 0.9 mg/ml. These peculiarities 
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Table I. Fit parameters <f>, \Z\, and A/Aq, for the BSA con- 
centration series without added NaCl, with intensities shown 
in Fig. 2. The additional parameters Lb — 0.711 nm, 
a — gb 2 = 7.40 nm, a = 1.75 nm, b = 4.74 nm are kept 
fixed, and Aq is taken from the form factor fit in Fig. 1. The 
obtained fit values at c p = 0.9, 1.8, 45, and 90 mg/ml should 
be taken with a pinch of salt (see related text). The listed 
values for c p are according to Eq. (2). 



c p [mg/ml] <f> |Z| n, [fiM] A/A 



0.9 


5.19 x 10~ 4 


34.5 


1216 


1.20 


1.8 


1.34 x 10" 3 


18.8 


608 


1.08 


4.5 


3.72 x 10~ 3 


19.1 


1278 


0.96 


7.2 


6.97 x 10" 3 


16.7 


1497 


0.97 


9 


1.04 x 10" 2 


14.6 


1510 


1.05 


13.5 


1.28 x 10 -2 


12.6 


1297 


0.81 


18 


2.06 x 10" 2 


10.8 


1292 


0.85 


45 


8.19 x 10~ 2 


9.47 


2375 


1.0 


90 


1.74 x 10 _1 


8.52 


3323 


1.0 



can be attributed to impurity contributions neglected in 
Eq. (8). Note also that the maximal intensities in both 
systems occur at wavenumbers well below the value 0.3 
nm -1 , where impurities are found to obstruct also the 
form factor fit in Fig. 1. 

All our attempts to remedy these fitting problems for 
the two most dilute samples failed. Lacking information 
about the shape and size distribution, and the interac- 
tions of the impurities, we cannot improve on Eq. (8). Re- 
stricting the wavenumber interval in the fitting procedure 
to q > 0.3 nm -1 leads to no improvement, either. While 
Eq. (8) is expected to be quite accurate in this restricted 
g-range, the maximum in I(q) is not included. The in- 
tensity for q > 0.3 nm -1 is a monotonically decaying 
curve, almost completely determined by the form factor. 
It therefore lacks distinct features coming from particle 
correlations, rendering the fit with respect to {\Z\, n s , (j)} 
into an overdetcrmined problem. For all these reasons, 
our fit parameters in Table I for c p = 0.9 and 1.8 mg/ml 
should not be considered as quantitatively accurate. 

Except for the two most dilute systems, all other sys- 
tems with concentrations from c p = 4.5 to 90 mg/ml in- 
cluded in Fig. 2 can be excellently fitted by Eq. (8). The 
obtained effective charges, salt concentrations, and vol- 
ume fractions all assume reasonable values, showing sys- 
tematic dependencies on the BSA concentration. Note, 
however, that for c p = 45 and 90 mg/ml, the SY model is 
pushed to its limit. On assuming a Hamakcr constant of 3 
ksT [67], the repulsive barrier height of the DLVO poten- 
tial becomes very small, with values of 1.3 and 0.5 ksT at 
c p = 45 and 90 mg/ml, respectively. The contact value of 
g(x) at x = 1 is just barely zero for the more dilute sys- 
tem, whereas g{x = 1 + ) w 0.9 in the more concentrated 



system. Obviously, the SY model with purely repulsive, 
spherically symmetric pair interactions is bound to fail 
when the particles are allowed to come into hard-core 
contact. Thus, the system with c p = 45 mg/ml, and fit- 
ted volume fraction <j> = 8.19%, is clearly on the border- 
line of the SY model. Somewhat unexpectedly, and prob- 
ably fortuitously, the system with c p = 90 mg/ml can still 
be fitted with good accuracy. Summarizing, the fit val- 
ues for the most concentrated systems with c p = 45 and 
90 mg/ml in Table I should be interpreted with caution, 
since the fit parameters might be significantly distorted 
by the discussed deficiencies of the SY model. An indi- 
cation for this could be the obtained fit values for 4>(c p ), 
which for the two most concentrated samples clearly over- 
shoot the linear dependence on c p found approximately 
for the lesser concentrated systems (see Table I). 

In closing our discussion of the static scattered inten- 
sities, we note that fit parameters slightly different from 
the ones in Table I are obtained, when in place of the BSA 
model spheroid axes (a, b) = (1.75 nm, 4.74 nm), the val- 
ues (a, b) = (1.80 nm, 4.60 nm) given in [9] are used. For 
instance, at c p = 4.5 and 18 mg/ml, the best-fit values 
for |Z| change to 18.4 and 10.7, respectively. Note that, 
in comparison to [32] , where the RMSA was employed in 
fitting I(q), we use here the improved MPB-RMSA inte- 
gral equation scheme for S(q), resulting in more precise 
fit- values. Moreover, different from the earlier intensity 
fitting described in [32], the dephasing influence on I(q) 
originating from the particle asphericity is accounted for 
approximately in the decoupling approximation used in 
the present study. The slightly different spheroid semi- 
axes (a, b) = (1.80 nm,4.60 nm), and the corresponding, 
slightly changed fit-parameters, do not cause appreciable 
changes in the dynamical properties. For instance, the 
collective diffusion coefficient changes by no more than 
3%, and the changes in the static- and high-frequency 
viscosities are less than 0.1%. Note that the somewhat 
smaller spheroid causes changes of the fitted volume frac- 
tion of about 5% which does not change absolute values 
but slightly rescales the protein concentration axis for the 
theoretical predictions. 



V. DYNAMIC PROPERTIES: EXPERIMENT 
AND THEORY 

In the following, we compare the DLS data for the 
collective diffusion coefficient of BSA solutions, and the 
static shear viscosity measured in our suspended couctte- 
type rheometer, to the results of the dynamic schemes 
discussed in Sec. III. Moreover, we test the validity 
of a generalized Stokes-Einstein relation connecting the 
viscosity to the collective diffusion coefficient and the 
isothermal osmotic compressibility. We reemphasize here 
that the employed theoretical schemes use S(q) and g(r) 
as the only input. With S(q) and g(r) determined from 
the fits to the SAXS-intensities, all theoretical results for 
dc, Voo and r\ are thus obtained without any additional 



11 



adjustable parameters. 



Collective diffusion coefficient 



Fig. 



(upper part) and d c 



3 includes our SLS/DLS data for l/S{q -> 0) 
dc (lower part), for aqueous 
BSA solutions in comparison with the theoretical pre- 
dictions. Systems without added salt, and for concentra- 
tions n s = 5 and 150 mM of added NaCl, are considered. 
Additional measurements using 500 mM of added NaCl 
(data not shown) agree almost perfectly with the data for 
n a = 150 mM, indicating that electrostatic repulsion is 
fully screened already at n s = 150 mM. As the input to 
the dynamics schemes, S(q) and g(r) were generated by 
the MPB-RMSA, using concentration-interpolated input 
parameters 4>(c p ) and Z(c p ) based on Table I. For no 
added salt, n s (c p ) was interpolated using Table I, while 
n a = 5 and 150 mM were kept fixed (independent of c p ) 
in the corresponding theoretical calculations. The value 

cIq 1 = 5.82 A /ns of the spheroid translational free diffu- 
sion coefficient was used to obtain dc in the experimental 
units from the dimensionless results for dc/do obtained 
by both theoretical schemes. 

For no added salt, the experimental dc{c p ) assumes a 
maximum at c p m 10 — 20 mg/ml. This maximum is qual- 
itatively reproduced by both theoretical schemes (cor- 
rected (57 and PA), but its location is predicted to occur 
at somewhat larger concentrations c p ~ 20 — 30 mg/ml. 
For BSA concentrations larger than the concentration at 
the maximum value for dc, the PA-predicted dc{c p ) re- 
duces strongly, eventually reaching unphysical negative 
values for c p > 110 mg/ml. This illustrates the expected 
failure of the PA scheme at higher concentrations, indi- 
cating that three-body contributions to HI, totally left 
out in the PA, but not in the 8^/ scheme, come into play 
f° r c p ^ 30 mg/ml. Up to the concentration value at 
the maximum of dc, both schemes agree very well, with 
residual differences not visible for c p < 20 mg/ml on the 
scale of Fig. 3. Despite its residual small inaccuracies, the 
self-part corrected expansion will therefore be used in 
the following calculations of dc- 

The physical origin of the non-monotonous concentra- 
tion dependence dc{c p ) at low concentrations of salt can 
be understood on the basis of Eq. (15), rewritten using 
dc 



d c as 



dc ,. H(q) 
— = hm — . 

do 9->o b(q) 



(20) 



with H(q) = d s /d + H d (q). The ratio in Eq. (20) con- 
sists of two competing factors. The factor l/S(q — > 0), 
inversely proportional to the isothermal osmotic com- 
pressibility of ideally monodisperse particles, increases 
monotonically as a function of the BSA concentration. 
Owing to the larger coupling constant 7 in Eq. (12a), a 
much steeper initial increase of \/S{q — > 0) is observed 
for weakly screened systems than for systems with added 
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Fig. 3. Top panel: Inverse zero-wavenumber limiting static 
structure factor of BSA solutions, obtained from SLS (con- 
nected black circles) and our MPB-RMSA scheme (red solid 
lines). Number concentrations, n s , of added NaCl as indi- 
cated. Bottom panel: Fast mode coefficient, Di = d^., ob- 
tained from the discussed double-exponential fit to the DLS 
data of BSA solutions (connected black circles), and d§ cal- 
culated by the self-part corrected 57 scheme (red solid lines), 
and the PA approximation (blue dotted curves). All theoret- 
ical curves are based on input parameters <f>(c p ) and Z(c p ) 
interpolated from Table I. In the zero added-salt case, the 
n„(c p ) values were also interpolated using Table I. Theo- 
retical results for added NaCl are obtained using fixed salt 
concentrations of n s = 5 and 150 mM. The input parameters 
L B = 0.711 nm, a = a B2 = 7.40 nm, a = 1.75 nm, b = 4.74 
nm, and do = df(a,b) = 5.82 A / ns are kept fixed through- 
out. For the zero added-salt case, the green vertical line seg- 
ment at c p « 34 mg/ml marks the protein mass-concentration 
where the surface-released counterion contribution to k 2 in 
Eq. (12b) is equal to the coion contribution. 



salt (c.f. the top panel of Fig. 3). As c p is further in- 
creased, the amount of surface-released counterions in- 
creases correspondingly, leading to an enhanced electro- 
static screening. As a consequence, the rate of change 
of l/S(q —> 0) with c p reduces significantly at a colloid 
concentration roughly set by the criterion, k^(c p ) = kg, 
of equal surface released counterion and salt-co-ion con- 
tributions to the screening parameter in Eq. (12b). 

The nominator in Eq. (20) is the reduced sedimenta- 
tion velocity, H{q — > 0), which is known from theory 
and experiment [53] to decrease monotonically, for not 
too large concentrations and low salinity according to 
1 — a sed 1 ^ 3 , w hh a sec i = 1.6 — 1.8 in the case of highly 
charged particles, and as 1 - 6.546 + 21.918 2 + 0(</> 3 ) 
for neutral hard spheres [68]. For strongly correlated 
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particles, the competition between decreasing compress- 
ibility and decreasing sedimentation coefficient with in- 
creasing c p leads thus to a maximum in dc{c p ), at a con- 
centration roughly determined from fc^(c p ) = k^. 

The DLS-measured values for dc are not quantitatively 
reproduced by the self-part corrected £7 scheme. Both 
in the zero added-salt case, and for n s = 150 mM, dc 
is underestimated by the corrected Sj scheme prediction 
by about 25%. The difference might be simply due to 
the complex-shaped BSA proteins having a translational 
free diffusion coefficient larger than the value dg" =5.82 

A 2 /ns used in the SY model. In fact, an extrapolation of 
the experimental data for dc to zero concentration leads 

to a larger value for do in the range of 6 — 7 A /ns, which 
can completely explain the differences in dc between ex- 
periment and theory. However, this low-concentration 
extrapolation should not be over-interpreted as being 
conclusive, since the experimental data are rather noisy 
for low concentrations. 

While the agreement between the theoretical and the 
experimental dc's is overall rather satisfying for very low 
and very high salt content, strong differences are found 
for the intermediate added NaCl concentration of 5 mM. 
This is not surprising, however, since already the zero 
added-salt experiments led to fit values for n s of 1 to 3 
mM. Therefore, n s is most probably a function of c p also 
in the 5 mM added NaCl case, instead of being constant 
as assumed in the calculations. Moreover, there is no 
obvious reason to expect that the relation Z(c p ), inter- 
polated from Table I, remains valid at arbitrary added 
salt concentrations. Additional future SAXS measure- 
ments at 5 mM added NaCl are necessary to determine, 
for this case, the precise dependence of n s and Z on c p . 

B. Static viscosity 

The rheometric results for 77 without added salt, and 
with 150 mM of added NaCl, are plotted in Fig. 4 as a 
function of c p , and compared to the theoretical predic- 
tions. Apart from pronounced differences at lower con- 
centrations, discussed in detail further down, the experi- 
mental data agree overall decently well with the theoret- 
ical predictions. Due to the rather weak microstructural 
ordering of the BSA proteins, characterized by structure 
factor peak heights less than 1.2 even for the most con- 
centrated samples, the shear-stress relaxation term A77 
contributes only little to 77, with a maximum relative con- 
tribution of about 10% near c p = 100 mg/ml. The domi- 
nant contribution to r\ is given by r?^ , which is predicted 
to good accuracy both by the PA scheme and the cor- 
rected £7 scheme, with practically equal results. The PA 
scheme is applicable to the whole experimentally probed 
concentration range of c p < 100 mg/ml, since three-body 
and higher order HI contributions affect rjoo to a lesser 
extent than dc (c.f. here Fig. 3, showing the failure of 
the PA prediction for dc already for c p < 50 g/l). 

The addition of larger amounts of salt lowers the val- 
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Fig. 4. Static relative (top panel) and reduced (bottom 
panel) viscosity for T = 25° C as function of c p . Theoret- 
ical curves are based on input parameters <j>(c p ) and Z(c p ), 
concentration-interpolated using Table I. Symbols: exper- 
imental data without added salt (black circles) and with 
n s = 150 mM (red diamonds). Lines: theoretical results 
without added salt (black solid line, n s (c p ) interpolated using 
Table I) and with a fixed salt concentration of n s = 150 mM 
(red dashed curve). Note the different c p -ranges in the two 
panels of the figure. 



ues for 77, as can be noticed from the two experimental 
data sets depicted in Fig. 4. The reason for this is the 
enhanced electrostatic screening, causing A77 to decrease 
with increasing salinity in going from strongly structured, 
charged spheres to basically neutral hard spheres. In con- 
trast, rjoa is known from theory and experiment [17] to 
increase upon the addition of salt, due to the enlarged 
influence of near-field His when the particles are allowed 
to get closer to each other in electrostatically screened 
systems. Thus, 77^ and A77 have opposite trends in their 
dependencies on the concentration of added salt. These 
competing trends are the reason for the weak crossover in 
the two theoretical curves for 77, noticed in the top panel 
of Fig. 4 at c p w 67 mg/ml. For particle concentrations 
larger than this c p value, the increase of 77^ overcom- 
pensates the decrease in A77 when, in place of the zero 
added-salt system, a system with n s = 150 mM is con- 
sidered. That such a weak crossover is not observed in 
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the experiment data in Fig. 4, points to an underesti- 
mation of the crossover concentration by our simplifying 
theories for r), possibly due to the neglect of His in the 
Arj calculation. 

A remarkable feature is noticed from the bottom panel 
of Fig. 4, where we plot the so-called reduced viscosity, 



as a function of c p . The function c p r] re< i/(j> reduces to the 
intrinsic viscosity, [77] at very low volume fractions where 
V ~* Vo + [v\4>- Features of dilute systems are more clearly 
revealed in r) re d than in 77. 

Both experimental data sets in the bottom panel of 
Fig. 4 show a local maximum of r\ re d at low c p values, 
which for the zero added-salt system (black open circles) 
is visible as a weak non-monotonicity near c p ~ 3 mg/ml. 
For the system with 150 mM added NaCl (red open di- 
amonds), the experimental maximum is represented es- 
sentially by a single data point at c p = 1 mg/ml, where 
Vred ~ 6.5 ml/g, whereas the remaining data points de- 
scribe a nearly constant plateau value of 4.5 ml/g. This 
plateau value is in good overall agreement with reported 
values for rj re d at low c p , in the range of 3.8 to 4.9 ml/g 
[69-72]. 

Regarding the large experiment error bars at very low 
c p , from the figure, we can not attribute physical signifi- 
cance to the single-point maximum in the n s = 150 mM 
system. A more refined data resolution in a future exper- 
imental study is clearly needed here. Even the maximum 
in r) re d for the zero added-salt case might be disputable 
on basis of the experimental data alone. However, the 
existence of such a maximum in r) re d draws its credibility 
from the comparison to the theoretical results, showing 
a maximum in r] re d(c p ) at a slightly lower value of c p . A 
similar non-monotonic behavior of rj re£ z(cp), with a pro- 
nounced peak at low c p , has been measured also in poly- 
electrolyte systems [73-75] , in low-salinity suspensions of 
charged silica spheres [76], and in microgels [77]. The 
effect has been described theoretically by scaling argu- 
ments [78], by the Rice-Kirkwood equation [79] for the 
shear viscosity in combination with a screened Coulomb 
potential [80], and for rod- like particles using a MCT 
scheme similar to ours [81]. In these earlier treatments, 
HI has been disregarded altogether. In our approach, HI 
is included in the for the present systems dominating rj^ 
part of rj. 

To rule out that the non-monotonicity of the theoret- 
ical r]red(c p ) is caused by BSA-specific dependencies of 
\Z\ and n s on c p , (c.f. Table I), we have investigated ad- 
ditionally a model system for fixed \Z\ = 20 and n s = 1 
mM, where we find again a maximum in r) re d(c p ). Thus, 
the maximum in rj re d{c p ) is a generic effect in weakly 
screened HSY fluids. It is entirely due to the shear-stress 
relaxation term A77, for (77^ — 770) / (^o^O increases mono- 
tonically in c p at arbitrary salt concentration. Since the 
His are neglected in our MCT treatment of the shear- 
stress relaxation part A77, we conclude that the local 



maximum in rj re d is basically a non-hydrodynamic effect, 
arising from electrostatic repulsion. We point out that 
the discussed physical mechanism underlying the non- 
monotonic behavior of rj re d{c p ) is different from the one 
causing the maximum in dc as a function of c p . The lat- 
ter maximum originates from a competition between elec- 
trostatic repulsion and hydrodynamic slowing in crowded 
systems. It is therefore not surprising that the maxima 
in rj re d and dc are located at considerably different pro- 
tein concentrations. Whereas the maximum of dc occurs 
at c p ~ 30 mg/ml (c.f. Fig. 3), the maximum in r\ re d is 
observed at c p < 5 mg/ml. 

The theoretical values for r\ re d in Fig. 4 underesti- 
mate the experimental data by a factor of about 1/2. 
In the low-concentration regime, the theoretical result 
for i] re dC p /4> approaches [77] = 2.5, owing to the under- 
lying effective sphere model. The intrinsic viscosity of 
BSA modeled as a spheroid is [rff 11 = 3.25, which is 
larger than the value for a sphere by a factor of 1.3 only. 
Therefore, this can not be the only cause for the observed 
deviation. However, the actual intrinsic viscosity of a 
heart-like shaped BSA protein is neither equal to that of 
a spheroid nor to that of an effective sphere. We recall 
here our discussion of Fig. 3, where we argued that do for 
a BSA protein might well be about 25% larger than the 
free diffusion coefficient, e?§", of the model spheroid. We 
can similarly argue that the observed differences between 
the experimental and theoretical r\ re d may be largely due 
to a value for the intrinsic viscosity of BSA of about 4 — 5, 
which is 20 — 50% larger than [r]] eU , and about twice as 
large as the [77] value of spheres. This could explain the 
observed difference. 

Finally, we note here that electrokinetic contributions 
to 77, ds and dc, originating from the non-instantaneous 
response of the microion- clouds around each protein, are 
not included in our treatment. Microion electrokinet- 
ics has the effect of lowering somewhat the values of ds 
and dc [61, 82], while enlarging the viscosity 77 [83, 84]. 
These effects can be expected to be stronger when kT 1 
is approximately equal to the particle size. Furthermore, 
electrokinetic effects are expected to be less significant at 
higher protein concentrations [85, 86]. 

C. Relation between viscosity and collective 
diffusion 

Kholodcnko and Douglas [19] have proposed the ap- 
proximate generalized Stokcs-Einstcin (GSE) relation 

(22) 

between collective diffusion coefficient, (static) viscosity 
and the square-root of the isothermal osmotic compress- 
ibility coefficient S(q — > 0,0). If this relation were ex- 
actly valid, the dimcnsionless function on the left-hand- 
side (lhs) of Eq. (22) would be a constant equal to one. 
The (approximate) validity of a GSE relation would be 
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c p [mg/ml] 

Fig. 5. Test of the long-time and short-time KD-GSE rela- 
tions in Eq. (22), with rf — rj and 77* = rjoo, respectively. Re- 
sults for BSA solutions without added salt (upper datasets), 
and with 150 mM of added NaCl (lower datasets) are shown. 
Red symbols: combination of d c from DLS, 77/770 from sus- 
pended couette rheometry, and S(q — s> 0) from SLS. Black 
lines: Theoretical results, combining do ~ d c and 7700 calcu- 
lated from the self-part corrected £7 scheme with S(q —¥ 0) 
from the MPB-RMSA scheme. For the long-time GSE ver- 
sion, 77 = r/oo + A77, with A77 from MCT is used. Lower bound- 
aries of the theoretical curves correspond to the short-time 
GSE, upper boundaries to the long-time version. The theo- 
retical curves are based on 5'(q , )-input with 0(c p ) and Z(c p ) 
concentration-interpolated using Table I. For zero added salt, 
n 3 (cp) was also interpolated on basis of Table I. The param- 
eters Lb = 0.711 nm and a = <tb 2 = 7.40 nm are kept fixed. 



very useful from an experimental viewpoint, since it al- 
lows to infer viscoelastic properties such as -q^ and 77 
from a dynamic scattering experiment where diffusion 
coefficients are determined. This is of particular rele- 
vance when the amount of protein available is too small 
for a mechanical rheometer measurement. Since we have 
experimental data sets for 77, dc, and S(q — > 0) for BSA 
solutions with low and high salt content at our disposal, 
together with theoretical tools to calculate these proper- 
ties, we are in the position to scrutinize the accuracy of 
the KD-GSE relation. We can do this not only for the 
special case of BSA solutions, but with our theoretical 
methods more generally for arbitrary spherical colloidal 
particles interacting by the HSY potential in Eq. (11). 

In their discussion of the GSE relation in Eq. (22), 
based on mode-coupling theory like arguments, Kholo- 
denko and Douglas have considered explicitly a dilute 
suspension of colloidal hard spheres to first order in <p 
only, where 7700 and 77 arc identical, since A77 = O(0 2 ). 
For high concentrations, we test now the validity of both 
the long-time and short-time versions of the KD-GSE 
relation, on recalling that different from 7700 and 77, c? c 
and cLq are practically equal even at high concentra- 
tions. In Ref. [19], it was argued that for uncharged 
hard spheres (HS) the KD-GSE relation is valid to lin- 
ear order in <p. We can check this statement analytically 



using numerically precise 2 nd order virial expansion re- 
sults for d% s = (d,S) HS , tj^, 5 , 7 ? ffS [68, 87, 88], and witl 
S HS (q — > 0,<f>) calculated from the precise Carnahan 
Starling equation of state. In this way, we obtain 



^Js HS {q 0) = 1 - 0.0460 + 1.37130 2 + O(0 3 ), 

(23a) 

^— '- — \/S HS (q->0) = 1 -0.0460 + 2.2820 2 + O((/> 3 ). 

(23b) 



dovo 



The short- and long-time versions of the KD-GSE rela- 
tion for hard spheres are identical to linear order in <p, 
with a coefficient, —0.046, which is not precisely vanish- 
ing but close to zero. However, to quadratic order in </> al- 
ready, where particle correlations come into play and 7700 
needs to be distinguished from 77, both GSE variants have 
distinctly non-zero virial coefficients. Since precise val- 
ues for the higher-order virial coefficients are not known 
to date, a test of Eq. (22) for larger <p can be made only 
using simulation and experimental data for dc(4>), ?7oo(</0 
and r](4>). This test has been performed in [8], where it 
is shown that both variants of the KD-GSE relation are 
approximately valid for hard spheres for </> < 0.1 only. 

Since neutral hard spheres are a special case of the 
HSY model, attained for 7 = or fc — > oo, as a matter 
of principle the validity of the KD-GSE relation for HSY 
systems is disproved already at this point. However, it 
still remains to be investigated in which concentration 
range the two KD-GSE relations are significantly violated 
when, instead of neutral hard spheres, weakly screened, 
charged HSY-likc particles such as charged proteins are 
considered. Note here that a virial expansion cannot be 
reasonably applied to charged particles at lower salinity, 
since the pair structure functions and thermodynamic 
properties in these systems depend on </3, 7 and A: in a 
non-analytical way. 

In Fig. 5, we plot the lhs function in Eq. (22), both 
in its short- and long-time form, as a function of c p . 
Both BSA solutions without added salt, and solutions 
with n s = 150 mM arc considered. Apart from a con- 
stant factor, which is related to the actual value of e?o 
in BSA solutions discussed earlier, the theoretical curves 
compare reasonably well to the experimental data. There 
are only small differences in the short-time and long-time 
GSE curves in the case of BSA solutions. 

With the hard-sphere-like behavior of the particles 
practically reached for n s = 150 mM, in the added-salt 
system the two KD-GSE variants apply for concentra- 
tions up to c p ~ 50 mg/ml, corresponding to « 0.1. For 
more concentrated systems, the lhs function in Eq. (22) 
increases initially, going trough a shallow maximum near 
c p « 90 mg/ml. For zero added salt, violation of the KD- 
GSE relations is observed theoretically at all non-zero 
concentrations, and can be noticed in our experiment al- 
ready for c p < 1 mg/ml. 

In our discussion of the KD-GSE relation, we proceed 
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Fig. 6. Test of the short- and long-time versions of the KD- 
GSE relation in Eq. (22) for volume fractions from very dilute 
values to <f> = 30%, and various salt concentrations as indi- 
cated. The collective diffusion coefficient, d^/do ~ d^/do, 
and high-frequency limiting viscosity, rjoo/rjo, are obtained 
from the self-part corrected scheme. The static viscos- 
ity, r\ — rjao + Ar/, is calculated using MCT for A77. Values for 
S{q -¥ 0) are obtained from the MPB-RMSA method. Input 
parameters Lb = 0.711 nm, a — 7.40 nm, and \Z\ = 10 are 
kept constant. 



now by characterizing the crossover behavior in going 
from the low-salt to the high-salt regime. To this end, in 
Fig. 6, we plot the lhs of Eq. (22) as a function of <j> for 
various salt contents, using the parameters Lb = 0.711 
nm, a = 7 AO nm, and \Z\ = 10. These parameters 
are typical of aqueous solutions of small globular pro- 
teins such as BSA, Lysozyme [11] and Apoferritin [61]. 
The charge number Z is kept constant here for simplic- 
ity. Theoretical results are plotted as a function of (f> 
instead of c p . In lowering the salt content in Fig. 6 step- 
wise by factors of 0.1, starting from a maximal value of 
n s = 100 mM, we find that the maximal (positive) devi- 
ation from one of the lhs function in Eq. (22) increases 
roughly logarithmically. For low salt content, n s < 1 
mM, the physical origin of the maxima in Fig. 6 is un- 
derstood from comparing the theoretical results for dc 
and r\ in Figs. 3 and 4, respectively: The maximal viola- 
tion of the KD-GSE relations occurs roughly at a volume 
fraction where dc(<t>) attains its maximum, i.e. for <fi de- 
termined approximately from k% ((f)) = k 2 s . Recalling that 
k\ oc (p and k 2 oc n s , this explains why the (fy location 
of the maxima in Fig. 6 shows a power-law dependence 
on n s for n s < 1 mM. For larger n s , a crossover to hard- 
sphcrc-likc behavior occurs, where the KD-GSE relations 
apply for <fi < 0.1. 



VI. CONCLUSIONS 

We have investigated static and dynamic properties 
of aqueous BSA solutions in an integrated conceptual 



framework, combining SLS/DLS, SAXS, and rheomet- 
ric measurements with analytical colloid theory. Solu- 
tions with physiological concentrations of added NaCl 
have been studied, as well as low-salt solutions showing 
distinct features in the concentration-dependence of the 
collective diffusion coefficient and the (reduced) viscos- 
ity. In our analytical theoretical approach, we have used 
a simple spheroid- Yukawa model of BSA with isotropic, 
repulsive pair interactions to calculate the static scat- 
tered intensity using the efficient MPB-RMSA method 
in combination with the orientational-translational de- 
coupling approximation. The form factor fit has been 
kept intentionally simple, without expecting extreme ac- 
curacy. The resulting S(q) have been used without any 
further fitting, in calculating dc, ?7oo, and 77 on basis 
of our well-tested theoretical methods. We have used 
the spheroid- Yukawa model for I(q), and the related ef- 
fective sphere- Yukawa model for the dynamic properties 
as minimal models without including additional protein- 
specific features, to clearly reveal the pros and cons of 
the model. This should help to point out more clearly 
the significance of left-out protein specific features. 

The measured static and dynamic properties of BSA 
are captured reasonably well in our simplifying SY model, 
with at least semi-quantitative accuracy, for mass concen- 
trations up to c p « 100 mg/ml. In the range 2 mg/ml 
^ c p < 50 mg/ml, reliable values for the effective protein 
charge number, and the residual electrolyte concentra- 
tion, have been obtained from the fits to the SAXS in- 
tensities. The SAXS fits are considerably obstructed for 
c p < 2 mg/ml by the presence of scattering impurities, 
and by the breakdown of the decoupling approximation 
f° r c p 50 mg/ml. 

A well-developed maximum in the concentration de- 
pendence of the collective diffusion coefficient of BSA 
was found at low salinity. This behavior is seen also 
in charge-stabilized colloidal suspensions. It is caused 
by the competition between electrostatic repulsion and 
hydrodynamic slowing down in crowded systems. More- 
over, a non-monotonic concentration dependence of the 
reduced viscosity of low-salinity BSA solutions was pre- 
dicted theoretically, and to some extent also seen ex- 
perimentally We have explained the local maximum in 
Vred(c p ) as a basically non- hydrodynamic effect caused 
by electric repulsion. A non-monotonic concentration- 
dependence of r] re d, with a pronounced peak at low con- 
centration, is observed also in polyelectrolyte solutions. 
Thus, the low-c p peak in r\ re d is a generic feature of 
charge-stabilized dispersions at low salinity. 

An essentially concentration-independent underesti- 
mation of the experimental dc and r\ re d by about 25% 
and 50%, respectively, is made in the theoretical predic- 
tions. Possible reasons for this are impurity effects, and 
an underestimation of the corresponding single-particle 
coefficients do and [77] through our disregarding of the 
complex protein shape and hydration shell morphology. 

We have analyzed the validity of a GSE relation by 
Kholodcnko and Douglas [19], which connects the col- 



16 



lectivc diffusion coefficient to the shear viscosity and to 
the isothermal osmotic compressibility. Despite its ap- 
pealing simplicity, the KD-GSE relation fails to capture 
the essential richness of macromolecular collective diffu- 
sion. It applies to decent accuracy to electrostatically 
screened solutions at high salinity, for volume fractions 
up to about 0.1. However, it is violated for more crowded 
high-salt solutions, and for basically all volume fractions 
under low-salt conditions. 

Possible extensions of the present work, which allow to 
maintain analytical simplicity to some extent, are the in- 
clusion of short-range attractive interactions for suspen- 
sions of larger salt content using, e.g. a two- Yukawa pair 
potential [89, 90], and the inclusion of surface patchiness 
[11]. For the static viscosity of more strongly concen- 
trated protein solutions than considered in the present 
work, the shear stress relaxation contribution, A?/, can 
become large in comparison to rj^. In calculating A?y, 
one needs then to account for HI contributions which 



tend to further enlarge its value. Such an inclusion of 
HI effects into A?y can be accomplished on basis of an 
extended MCT scheme discussed in Refs. [65, 66]. These 
extensions will be the subject of a future study. 
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